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Abstract. - We study the Casimir effect for the perfect Bose-gase in the slab geometry for 
various boundary conditions. We show that the grand canonical potential per unit area at the 
bulk critical chemical potential = has the standard asymptotic form with universal Casimir 
terms. 



In contrast to the well-known Casimir effect for the photon gas, see e.g. [1], this effect 
for the massive quantum particles [quantum gases) is much less explored. In the present 
letter we consider the case of the perfect Bose-gas. The perfect Bose-gas is the simplest 
quantum mechanical system showing the spontaneous breaking of a continuous symmetry 
(Bosc-Einstein condensation) for which the Casimir effect should appear. To our knowledge, 
this has surprisingly not been studied before. 

We first recall the definition of the Casimir amplitude at critical points or for phases with 
broken symmetry [2] [3]. Let 

^,(T)= hm ^$i(r,d) (1) 

be a limiting thermodynamic potential per unit area in the 3-dimensional slab geometry oo x 
oox d, where $ l (T, d) is the finite- volume thermodynamic potential in the box A = L x L x d 
and T denotes the temperature. The finite size scaling analysis in the critical regime shows 
that large d-asymptotics of ^ has usually the form: 

MT) = difbuikiT) + ^,ur f{T) + + ... . (2) 

Here ipbuik(T) is the potential density in thermodynamic limit and ipsurf{T) is the surface 
potential correction. The coefhcient A(T) in the third term in |j2Jl is called the Casimir 
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amplitude. The Casimir force between the slab faces is defined as 

F{d) = - ddiMT) - d<fbuik{T)] = 2 + . . . . (3) 

In this regime, i.e. at a critical point or in phases with a long-range correlations generated by 
the broken symmetry, the value of A(r) is expected to be non-zero and universal, depending 
only on the system and the boundary condition universality classes. Out of the critical regime 
the finite size corrections in slab are expected to be exponentially small, thus the Casimir 
amplitude A(r) = 0. 

The perfect Bose-gas grand- canonical thermodynamic potential in the box L x L x d has 
the form [4] 

ci>i,,(T,M) = ri^ln{l-e-^(^('')-'^)} , eik) ^^^(kl + kl + kl) , (4) 

k 

where P = (ksT)^^ and /i is the chemical potential. The sum in runs over the set defined 
by boundary conditions (b.c.) for the Laplacian operator in the box L x L x d. Below we 
shall consider (for simplicity) the periodic b.c. in the x-y directions: 

_ 27r _ 27r _ 

and three cases of b.c, Dirichlet [D), Neumann {N) and periodic (P), in the z direction: 

{D) fc^ = , = 1,2,... ; {N) fc, = ^n, , = 0, 1, 2, . . . (6) 

(P) fc, = ^n, , = 0,±1,±2,... . (7) 

Therefore, in the slab limit, L ^ oo, the grand-canonical thermodynamic potential per unit 
area is 

MT,,) - }^^^'^^AT,,) = j^l/\E^n{l-e-^l^^^^^^^^^ 

" ~ (2^ Jr^ Y e/'[^(q)+^(fc-)-fl - 1 ' *• 

where q ~ [q^, qy) is a two dimensional wave-vector in the (x, y) plane, £(q) — fi^iq^ + 9y)/2TO 
and e{kz) = h'^kl/2m. The last equality results from an integration by part with respect to 
9=|q|. 

It is well-known that in the bulk limit, L oo , d ^ oo, the Bose-Einstein condensation 
for the perfect Bose-gas may occur only for yu = 0. Therefore in the sequel we distinguish 
the two regimes : n < (normal phase) and /x = (condensed phase). Our results for the 
asymptotic expressions of ipd(T, /x) as d ^ oo are the following: 

(a) If /X < 0, then 

MT, /i) = -dpbulk(T, fi) + ^surf(T, fl)+0 (^e-^on^t-V^d/X^^ ^ (g) 

where 

Pbulk(T,^i) = J d^kln (l - e-^I^C^)-^!) (10) 
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is the standard Bosc-gas pressure and A = h^Jj3/rn denotes the thermal wave-length. We get 
that flulfC^^ l^) " ^^O'' periodic b.c., but 

Vsurf (T, M) = ±^ J d q-^^-^-p-^ (11) 

with + for Dirichlet and — for Neumann b.c. 



(b) If /i = 0; then we obtain for the periodic b.c. the asymptotics 



A(^)(r) 



-d Pbuik{T,Q) 
fcBTC(3) 



A(^)(T) + 0((A/d) 



for any M > 1 



(12) 



where Q{s) = X^i^i * is the Riemann zeta-function. For the Dirichlet and Neumann b.c. 
we find 



MT, 0) = -d ptuik{T,0) + ^[tfiT, 0) + ^ A(^'^)(T) + O {{X/df') 



fcBTC(3) 



, for any M > 1 
(13) 



and vi^;V{T, 0) is given by (EJ. 



We see that for Dirichlet and Neumann b.c. the Casimir amplitudes are the same, 
implying identical attractive forces in slab in presence of the Bose-Einstein condensation. 



To obtain these results we consider separately the cases ^ < and /i = 0. 
For /i < and Dirichlet b.c. (see ©) the result (a) follows from representation of as 
the sum: 

Since for < the function 4>{u) is infinitely differentiable at w = and all its derivatives 
vanish at infinity, we can use the Euler-MacLaurin theorem [5], [6] in the form 

oo poo 1 R H 

Y^F{n)= dsF{s)~-F{0)-^F^'\Q)~^F^-'\0)-..., (15) 

n=l 

where {-B2n}^i are the Bernoulli numbers and {F^'^^~^\s)}'^_^ denote the corresponding 
derivatives. When applying (|15|l to (|14|l . the two first terms of H15() yield respectively the 
bulk and surface contributions in Moreover, since 4>{u) — 4>{—u) all odd derivatives 

|0(2n-i) (y)|^^^ vanish at u = 0. So for large d the corrections to the first two terms are 
smaller than any power of c?~^. In fact, a more elaborated estimate gives an exponentially 
smah correction 0(e-c°"^* V^'^/^), see (HTJ below. 

If ^ = 0, all derivatives (/)(^"~^)(m) of order 2rt — 1 > 2 diverge at u = 0. So instead of the 
Euler-MacLaurin formula one has to use a more refined representation for the sum (|14|l . To 
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this end we use the Jacobi identity (see e.g. [7], Ch.ll): 



71=1 \ V /V ri—1 



which can be obtained as a simple consequence of the Poisson summation formula. First we 
represent (|14|l by its low activity series for /i < (setting v = /3e(q) = \^q^ /2) 



d«z;y y e'"-^e-'-["+(^^"/'^)'/2l = (17) 

^ "'U „=lr=l 



r—1 n—1 V V / r — 1 r—1 



2d A A 



where the last equation follows from the the Jacobi identity with a = rTT{X/d)'^ /2. These 
series are absolutely convergent when /i = and one recognises that the first term in the right- 
hand side of (|17|) in nothing but —d ptuikiT, 0) and the second one coincides with ^[^IfiT, 0). 
To obtain for // = the large d asymptotics of the last term we write the sum over r in the 
form {X/nd)^Sn{d)Sn{d), where 



f ^-2{nd/\f/r 



Sn{d) =J2^{6nid)r) , ^(rr) = ^ , S^id) - ( — ) . (18) 

r=l ^ 



Since the function tp^x) is infinitely differentiable for a; > and all its derivatives vanish at 
X = and at infinity, the Euler-MacLaurin formula (|15f) yields for large d and n 

/>oo / — 

5^{d)S^{d)^ dxi;ix)+0{{6,,{d))'')^^ + 0{{5^{d))'') for any A/ > 1 , (19) 

i.e. the corrections to the integral are smaller than any power of Sn{d). Inserting p9|) in the 
last term of the right-hand side of p7|) we get 



n=l 



, for any A/ > 1 



/3(V2^A)3^,^ r5/2 87r/3d2 

(20) 

that proves (|13|1 for the Dirichlct b.c. By © for Neumann b.c. sums in H16|) and 1)17(1 start 
with n = 0, which implies changing the sign of the surface term 1)11(1 and gives p3(l for 
Neumann b.c. 

For periodic b.c. Q we have to use instead of l(14() the representation 

oo oo 

MT,t^)= Hid/2r'n) = 2j2H{d/2)-'n) + m- 

n — — oo n—1 

By consequence, the contribution of the n = term, which corresponds to the surface correc- 
tions, will disappear and we obtain ((12(1 . 
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The Jacobi identity allows also to prove the asymptotic ©• To this end notice that 
double-sum in the last term of H17|) has for /i < the estimate 



r— 1 n— 1 r—1 n—1 

C(5/2) 



that proves Q for large d. 



0{e 



Now few remarks and comments are in order. 

The grand potential of a free Fermion gas does not have a Casimir term for any value of 
the chemical potential. Indeed replacing the Bose by the Fermi distribution in JSJ gives a 
corresponding (j>{u) function H14|l that is infinitely differentiable at u = for all ^. Since the 
function and all its derivatives vanish at u = the Eulcr-MacLaurin formula yields corrections 
smaller than any inverse power of d. 

According to common wisdom, the Casimir force is due to Goldstonc modes that will occur 
in the bulk limit when a continuous symmetry is spontaneously broken. Usually they are taken 
ad hoc as a part of phenomenological models, see e.g. [8] for the case of superfluid films. In 
the free Bose-gas it is explicit: the excitations e(k) — n become gapless when ^ — 0. Then by 
the London-Placzek formula [9], [10], this generates long-range particle-particle correlations 

p^'\r,,r2)-p' ^Po{T)\r,-r,\-\ (22) 

as |ri — oo in the condensed phase. Here po denotes the Bose-Einstein condensation 

density of the perfect Bose-gas. Casimir forces are usually attributed to such correlations in 
the critical regime. One can check that it is true for the perfect Bose-gas because of long-range 
correlations that appear in the slab limit. 



It is known that the perfect Bose gas particle number fluctuations in the condensed phase 
are abnormal, i. e. proportional to volume V in contrast to the normal thermodynamical 
fluctuations VV, because of lack of superstability, see e.g. [11]. So we would like to warn that 
the results for Casimir amplitudes might be very different from (I12II . H13|) in the presence of 
a superstable interaction. 

In the case of electromagnetic interactions, the Casimir term is always present as a result 
of the long range of the forces. In the standard calculation of the zero-temperature Casimir 
force between perfect conductors [1] the Casimir term appears because the third order deriva- 
tive occuring in the appropriate Euler-MacLaurin expansion does not vanish, contrary to the 
case (a) above. This is due to the linear form fkuk = fi-c|k| of the photon spectrum (non 
analytic at k = 0). Massive photons fiuk = cy/ (fi.|k|)2 + (mc)^ do not produce Casimir forces 
for the same reasons as for Fermions. On the other hand, in the high temperature regime, 
the Casimir term — C(3)/167r/3c?^ becomes purely classical and is entirely due to fluctuating 
Coulomb potentials in the conductor walls [12] - [14]. Moreover, a classical Coulomb system 
enclosed in a slab with perfectly conducting walls manifests an universal finite size correction 
of the Casimir type of the opposite sign C(3)/167r/3fi^ [15]. Such Coulomb systems are critical 
at all temperatures because of the potential fluctuations are always of a long range, although 
charge correlations are themselves of a short range as a consequence of screening. 
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Wc note that the Casimir terms found in (b) are classical and universal, namely they do not 
depend on the Planck constant and the particle mass. In fact, for the free Bose-gas, it follows 
from a simple dimensional analysis that such terms are necessarily of the form CkBT/cP 
where C is a numerical constant. They are present at all positive temperature provided that 
the density p of the gas is higher than the critical density Pc{P). Since M is arbitrary in l|12(l 
and H13() . we see that the Casimir amplitude is dominant as soon as d > A. For instance in 
liquid Helium at T = 2K, one has A - ^ 4A, i. e. a slab of thickness slightly large than 

the interatomic distance already shows the Casimir effect. If we fix d and let T ^ we enter 
in the regime A ^ d. By (|17fl it is sufficient to keep the ?7 = ?- = 1 term which gives 

So in contrast to the photon case the grand-potential is exponentially small as T — *■ 0, which 
excludes Casimir effect in this regime. 

* * * 

The paper was originated during V.A.Z.'s visits to Ecole Polytechnique Fcdcrale de Lau- 
sanne. He gratefully acknowledges warm hospitality at Institute of Theoretical Physics -EPFL. 
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